APPLICATIONS OF VARIATIONAL ANALYSIS 
TO A GENERALIZED HERON PROBLEM 



BORIS S. M0RDUKH0VICB3, NGUYEN MAU NA]V0 and JUAN SALINAS JR0 

Abstract. This paper is a continuation of our ongoing efforts to solve a number of geometric problems 
and their extensions by using advanced tools of variational analysis and generalized differentiation. Here we 
propose and study, from both qualitative and numerical viewpoints, the following optimal location problem 
, as well as its further extensions: on a given nonempty subset of a Banach space, find a point such that the 

sum of the distances from it to n given nonempty subsets of this space is minimal. This is a generalized 
version of the classical Heron problem: on a given straight line, find a point C such that the sum of the 
distances from C to the given points A and B is minimal. We show that the advanced variational techniques 
allow us to completely solve optimal location problems of this type in some important settings. 

Key words. Heron problem and its extensions, variational analysis and optimization, generalized differen- 
tiation, minimal time function, convex and nonconvex sets. 
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^ ■ 1 Introduction and Problem Formulation 

S3 . 

In this paper we propose and largely investigate various extensions of the Heron problem, which 
seem to be mathematically interesting and important for applications. In particular, the one of this 
type is to replace two given points in the classical Heron problem by finitely many nonempty closed 
\ subsets of a Banach space and to replace the straight line therein by another nonempty closed subset 

of this space. The reader are referred to our paper |14j for partial results concerning a convex version 
OO of this problem in the Euclidean space M n . 

| Recall that the classical Heron problem was posted by Heron from Alexandria (10-75 AS) in his 

Catroptica as follows: find a point on a straight line in the plane such that the sum of the distances 
from it to two given points is minimal; see [H [5] for more discussions. We formulate the distance 
function version of the generalized Heron problem as follows: 

n 

minimize D(x) := } d(x;Cli) subject to x £ O, (1-1) 

x, 

where fl and Qj, i = 1, . . . , n, n > 2, are given nonempty closed subsets of a Banach space X endowed 
with the norm || ■ ||, and where 

d(x;Q) :=inf{||a:-y|| | y G Q}. (1.2) 

is the usual distance from x £ X to a set Q. Observe that in this new formulation the generalized 
Heron problem (jl.ip is an extension of the generalized Fermat-Torricelli problem proposed and 
studied in [13]. The difference is that the latter problem in unconstrained, i.e., Vt = X in (|1.1|) 
while the presence of the geometric constraint in the generalized Heron version (|1.1[) makes it more 
mathematically complicated and more realistic for applications. Among the most natural areas of 
applications we mention constrained problems arising in location science, optimal networks, wireless 
communications, etc. We refer the reader to the corresponding discussions and results in |13j and the 
bibliographies therein concerning unconstrained Fcrmat-Torricelli-Steiner- Weber versions. Needless 
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to say that the presence of geometric (generally nonconvex) constraints in essentially changes 
these versions while referring us to the original Heron geometric problem. 

In fact, we are able to investigate a more general version of problem where the distance 

function (|1.2[) is replaced by the so-called minimal time function 



with the constant dynamics x 6 F C X and the target set Q C X in a Banach space X; see [12] 
and the references therein for more discussions and results on this class of functions important for 
various aspects of optimization theory and its numerous applications. 

The main problem under consideration in this paper, called below the generalized Heron problem, 
is formulated as follows: 



where F is a closed, bounded, and convex set containing the origin as an interior point, and where 
n and rti for i = 1, . . . , n are nonempty closed subsets of a Banach space X\ these are the standing 
assumptions of the paper. 

When F = IB in (|1.4[) . this problem reduces to the one in Note that involving the minimal 

time function (|1.3|) into (|1.4[) instead of the distance function in allows us to cover some 

important location models that cannot be encompassed by formalism cf. [T5] for the case of 

convex unconstrained problems of type (| 1 .4[) and |13j for the generalized Fermat-Torricelli problem 
corresponding to (|1.4|) with ft = X. 

A characteristic feature of the generalized Heron problem (|1.4|) and its distance function speci- 
fication is that they are intrinsically nonsmooth, since the functions (|1.2[) and (|1.3j) are non- 
differentiable. These problems are generally nonconvex while the convexity of both cost functions 
in and (|1.4|) follows from the convexity the sets Qj. This makes it natural to apply advanced 
methods and tools of variational analysis and generalized differentiation to study these problems. To 
proceed in this direction, we largely employ the recent results from [T^] on generalized differentiation 
of the minimal time function (jl.3l) in convex and nonconvex settings as well as comprehensive rules 
of generalized differential calculus. As can be seen from the solutions below, the constraint nature 
of the Heron problem and its extensions leads to new structural phenomena in comparison with the 
corresponding Fermat-Torricelli counterparts. Note that a number of the results obtained in this 
paper are new even for the unconstrained setting of the generalized Fermat-Torricelli problem. 

The rest of the paper is organized as follows. In Section 2, we present some basic constructions 
and properties from variational analysis that are widely used in the sequel. Section 3 concerns 
deriving necessary optimality conditions for solutions to the generalized Heron problem in the case of 
arbitrary closed sets and fi,, i = 1, . . . , n, in (|1.4[) and its specification The results obtained 

are expressed in terms of the limiting normal cone to closed sets in the sense of Mordukhovich 
[5]. We pay a special attention to the Hilbert space setting, which allows us to establish necessary 
(in some cases necessary and sufficient) optimality conditions in the most efficient forms. Some 
examples are given to illustrate applications of general results in particular situations. In Section 4 
we develop a numerical algorithm to solve some versions of the generalized Heron problem in finite 
dimensions while the concluding Section 5 is devoted to the implementation of this algorithm and 
its specifications in various settings of their own interest. 

Our notation is basically standard in the area of variational analysis and generalized differentia- 
tion; see [!|J[IB]. We recall some of them in the places they appear. 

2 Tools of Generalized Differentiation 

This section contains basic constructions and results of the generalized differentiation theory in 
variational analysis employed in what follows. The reader can find all the proofs, discussions, and 
additional material in the books [H [S3 [TOl HH1 [T7] and the references therein. 



T^{x) := inf {t > 0| Q n (x + tF) ^ 0} 



(1.3) 
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(1.4) 
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Given an cxtcndcd-real- valued function <p: X —> ]R := (— 00,00] with x from the domain 
dom<p := {x G X\ ip(x) < 00} and given e > 0, define first the e-subdifferential of tp at x by 

S.p(x) := {«• € JC hminf ^'^'f'"-^ > - e l (2.1) 

For e = the set d(p(x) := doip(x) is known as Frechet /regular subdifferential of <p at x. It follows 
from definition (|2.1I) that regular subgradicnts are described as follows: x* S d £ ip(x) if and only if 
for any r\ > there is 7 > such that 

(a;*, x — x) < tp(x) — tp(x) + (e + 77) ] I ac — x|| whenever .t G S + 7-2? 

with standing for the closed unit ball of the space in question. When ip is Frechet differen- 
tiable at x, its regular subdifferential d<p(x) reduces to the classical gradient {Vip(x)}. Despite the 
simple definition (|2.1I) closely related to the classical derivative, the regular subdifferential and its 
e-enlargcments in general do not happen to be appropriate for applications to the generalized Heron 
problem under consideration due to the serious lack of calculus rules. 

To get a better construction, we need to employ a certain robust limiting procedure, which lies 
at the heart of variational analysis. Recall that, given a set- valued mapping G: X X* between 
a Banach space X and its topological dual X*, the sequential Painleve-Kuratowski outer limit of G 
as x — > x is defined by 



LimsupG(x) := < x* G X* 



3 sequences Xk — > x, x* k ^> x* as k — > 00 

(2.2) 

such that x* k G G{x k ) for all k G N := {1, 2, . . .} j, 



where w* signifies the weak* topology of X*. Applying the limiting operation (|2.2[) to the set-valued 
mapping (x, e) =4 d £ (p(x) in (|2.1[) and using the notation x A x := x —> x with ip(x) — > tp{x) give us 
the subgradient set 

dip(x) := Limsup d e ip(x) (2-3) 

x — fx 

known as the M or -dukhovich /limiting subdifferential of ip at x. We can equivalently put e = in (|2.3[) 
if <^ is lower semicontinuous around x and if X is Asplund, i.e., each of its separable subspaces has a 
separable dual; the latter is automatics, e.g., when X is reflexive. Recall that tp is subdifferentially 
regular at x if dip(x) = dtp(x). 

Note that every convex function <p is subdifferentially regular at any point x G dom ip with the 
classical subdifferential representation 

dip(x) = {x* G X* I <x*,i- x) < ip(x) - <p(x) for all x G X}. (2.4) 

However, the latter property often fails in nonconvex setting, where dtp(x) may be empty (as for 
ip(x) = — \x\ at x = 0) with a poor calculus, while the limiting subdifferential (|2.3[) enjoys a /uZZ 
calculus (at least in Asplund spaces) due to variational/ extremal principles of variational analysis. 
We following calculus results are most useful in this paper. 

Theorem 2.1 (subdifferential sum rules). Let tp i ; X — > M, i = 1, . . . , n, be lower semicontinu- 
ous functions on a Banach space X . Suppose that all but one of them are locally Lipschitzian around 
x E (~\" =1 dom ipi . Then: 

(i) We have the inclusion 

n n 

(2.5) 

i-1 i=l 

provided that X is Asplund. Furthermore, inclusion (|2.5p becomes an equality if all the functions tpi 
are subdifferentially regular at x. 
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(ii) When all the functions (pi are convex, the equality 

n n 
i=l i=l 

ZioZcfe wit/i no Asplund space requirement. 

Note that assertion (ii) of Theorem I2.1[ which is the classical Moreau-Rockafellar theorem, is a 
consequence of assertion (i) in the case of Asplund spaces; see [91 Theorem 3.36]. 

Finally in this section, recall that the corresponding normal cones to a set Q at x £ f2 can be 
defined via the subdifferentials (|2.ip and (|2.3p of the indicator function by 

JV(x; ft) := 9<5(x; 0) and iV(sc; fi) := 3J(2;Q), (2.7) 
where <5(x; ft) = if x £ f2 and <5(x; Q) = oo otherwise. 



3 Optimality Conditions for the Generalized Heron Problem 

The main results of this section give necessary optimality conditions for the generalized Heron 
problem under consideration, which occur to be necessary and sufficient for optimality in the case of 
convex data. To begin with, we would like make sure that problem (|1.4[) admits an optimal solution 
under natural assumptions. 

Proposition 3.1 (existence of optimal solutions to the generalized Heron problem). The 

generalized Heron problem (|1.4[) admits an optimal solution in each of the following three cases: 

(i) X is a Banach space, and the constraint set Q is compact. 

(ii) X is finite- dimensional, and one of the sets fi and fli as i — 1, . . . ,n is bounded. 

(iii) X is reflexive, the sets Q and Cli as i = 1, . . . ,n are convex and one of them is bounded. 

Proof.lt follows from JTT1 Proposition 2.2] that the minimal time function (|1.3p and hence the 
function T in (|1.4[) are Lipschitz continuous. Thus the conclusion in the case (i) follows from the 
classical Weierstrass theorem. 
Consider the infimum value 

7 := inf T(x) < oo 

in problem (|1.4[) and take a minimizing sequence {xk} with T(xk) — > 7 as k — > 00 and Xk £ f2 for all 
k £ IN. Now assume that X is finite dimensional and Oi is bounded. When k is sufficiently large, 
one has 

T,i(x k ) <T(x k ) < 7 +l. 
Thus there exist < t k < 7 + 1, fk £ F, and Wk £ ^i such that 

Xk +t k fk = w k . 

Since both F and fii are bounded, (xk) is a bounded sequence, and hence it has subsequence that 
converges to x £ £1. Then a; is a solution of the problem under (ii). The proof in case (iii) is similar 
to that given in [T3J Proposition 4.1]. A 

To proceed with deriving optimality conditions for the generalized Heron problem (jl.4p and its 
specification (|1.1[) . we need more notation. Define the support level set 

C* := {x* £ X*\ a F (-x*) < 1} 

via the support function of the constant dynamics 

(Jf{x*) := sup(x*, re), x* £ X* . 
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The generalized projection to the target set Q via the minimal time function (|1.3|) is a set-valued 
mapping IIq : X X defined by 

HZ(x):=Qn{x + T$(x)F), x G X. (3.1) 
Considering further the Minkowski gauge 

p F {x) := inf {t > 0| x G tF}, x G X, (3.2) 
and involving the limiting normal cone from (|2.7[) . we define the sets 

r |J [-^(w-^nJV^Oi)] for ^fi^n^^^g, 

Mx) := I (3.3) 
I N(x; fii) n C* for igSI, as i = l,...,n. 

We say that the minimal time function Tq (•) is wdZ posed at a; if for every sequence {xfc} converging 
to x there is a sequence {wk} such that w^, G Hgixk) and {wk} contains a convergent subsequence. 
The reader is referred to [121 Proposition 6.2] for a number of verifiable conditions ensuring such a 
wcll-posedness of the minimal time function. 

Our first theorem establishes necessary as well as necessary and sufficient conditions for optimality 
in (|1.4[) via the sets Ai{x) from Q3.3P in general infinite-dimensional settings. 

Theorem 3.2 (necessary and sufficient optimality conditions for the generalized Heron 
problem in Banach and Asplund spaces). Given x G fi, suppose in the setting of (|1.4p that the 
minimal time function Tq, is well posed at x, for each i G {1, . . . , n] such that x ^ Oj. The following 
assertions hold: 

(i) Let x be a local optimal solution to ()1.4[) . and let X be Asplund. Then we have 

n 

eJ2Mx)+N(x;Q), (3.4) 

where the sets Ai(x) are defined in (|3.3[) . 

(ii) Let X be a general Banach space, and let all the sets fi and f2j as i = 1, . . . ,n be convex. 
Given x G 7 assume that TLq. (x) ^ for i = l,...,n with x £ f2j, select any uj G ff^. (x), and 
construct Ai{x) by 

Ai(x) := N(u>; Qi) n [ - dp F (ui - x)] for x <£ il, (3.5) 

and by the second formula in (|3.3p otherwise. Then x is an optimal solution to (|1.4[) if and only if 
inclusion (|3 .4[) is satisfied. 

Proof. Observe first that problem (|1.4p can be equivalently written in the form 

minimize T(x) + 8(x; CI). (3-6) 

It easily follows from definitions (|2.ip and (|2.3p of regular and limiting subgradients and their 
description (|2.4p for convex functions that the generalized Format rule 

G df(x) c df{x) (3.7) 

is a necessary condition for a local minimizer x of any function f: X —> 1R being also sufficient 
for this if / is convex. To justify now assertion (i), we apply (|3.7p via df(x) to the cost function 
f{x) := T{x) + Six; H) in (|3.6p and then use the subdifferential sum rule for limiting subgradients 
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from Theorem 12. lf i) in Asplund spaces by taking into account that the functions T£. are Lipschitz 
continuous. It follows in this way that 

ed(T + 5(--n))(x) c dT{x) + N{x-,n) 

n 

C^^(i) + ff(i;!l). (3 ' 8) 



<(5)C^(2), i = l,...,n. (3.9) 



Employing further the subdifferential formulas for the minimal time function from [13j Theorem 3.1 
and Theorem 3.2] gives us 

Sin 

Substituting the latter into (|3.8p justifies inclusion (|3.4[) in assertion (i) of the theorem. 

To justify assertion (ii), we apply Theorem I2.1f ii) for convex functions on Banach spaces and 
conclude in this way that both inclusions "C" in (|3.8[) hold as equalities and provide necessary and 
sufficient optimality conditions for optimality of x in (|1.4[) . Employing finally [T^l Theorem 7.1 and 
7.3] gives us the equalities in (|3.9[) . where the sets Ai{x) are calculated by (|3.5|) when x £ Oj. This 
completes the proof of the theorem. A 

It is not hard to check under our standing assumptions that the requirement 11^. (x) ^ in 
Theorem 13. 2f ii) is automatically satisfied when the space X is reflexive. 

The next theorem allows us to significantly simplify the calculation of the sets Ai(x) in Thco- 
rcm l3.2l for the case of Hilbert spaces and thus to ease the implementation of the optimality conditions 
obtained therein. Besides this, it leads us to an improvement of optimality under some additional 
assumptions. Namely, we can replace the limiting normal cone in (|3.4[) by the smaller regular one 
for an arbitrary closed constraint set fi. Define the index sets 

I(x) := {i G {l,...,n} | x G Oj} and J(x) = {i G {1, . . . , n] \ x g fy}, x G X. (3.10) 

We obviously have I{x) U J(x) = {1, . . . , n] and I(x) n J(x) = for all x G X. 

Theorem 3.3 (improved optimality conditions in Hilbert spaces). Consider version (|1.1[) 
of the generalized Heron problem with a Hilbert space X in the assumptions of Theorem 13.21 The 
following assertions hold: 

(i) Let x G f2 be a local optimal solution to (jl.lj) . and let IT(x; f2j) ^ whenever i £ J(x). Then 
for any Oi{x) G Ai(x) as i G J{x) we have 

- J2 a ^) e Mx)+N(x;n), (3.11) 

i£j(x) i£l{£) 

where each set Ai(x) is computed by 

,,- n \ for x £ 11,, 
Ai{x) = l a(x,Ui) (312) 

[ N(x;n t )nlB for x G fl, 
whenever i = 1, . . . , n. If in addition I(x) = 0, then 



J2 a i(%) GiV(x;Q). (3.13) 



i=l 



(ii) If all the sets fl and fli as i = 1, . . . , n are convex, then each set Ai(x) as i G J(s) Mi (|3.12|) 
is a singleton {ai(x)} and condition (|3.1ip is necessary and sufficient for the global optimality of 
x G O in problem (jl.ll) . 
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X — t/J' 

Proof. To justify assertion (i), pick Ui £ H(x;{li) for all i £ J(x) such that a%{x) — z and 

(tyX] \ t%) 

get the relationships 

n n 

Y \\x-Qi\\+ Yj d{x;n i )=Y,d(x;n i )<Y,d(x;n i )< Y ||a;-w i ||+ Y d(x;Sk) 

ieJ(x) iel(x) i=l i=l ieJ(x) iel(x) 

for all x € f2 around x. This shows that a; is a local optimal solution to the problem 

minimize p(x) := ||a; — a)j || + d(x;Qi) subject to x G f2. (3-14) 

i£j{x) i£l{x) 

Since the norm function on a Hilbert space is Frechet diffcrcntiable in any nonzero point, we conclude 
that each Pi(x) := \\x — Qi\\ as i £ is Frechet diffcrcntiable at a; with 

Vp l (5) = ip r^ji = -Tfi = Oi(x). 

— a(a;;S2ij 

Applying to (|3.14j) the first inclusion in the generalized Format rule (|3.7[) and then using the subd- 
iffcrential sum rules from [HI Proposition 1.107(i)] for regular subgradients and from Theorem 12. lf i) 
for limiting ones, we get 



Ed[p + S(-,n)](x) = J2 V Pt (x)+d[ Y d(-;Oi) + 5(-;fi)](s) 

i£j(x) i£l(x) 



eJ(x) ie-T(x) 

i£.J(x) i£l(x) 

c ^ Oi(x)+ ^ [N{x; Cli) nB] + JV(x; O) 
ie.j(x) iei(x) 

= Y a <(%)+ Y Mx)+N(x;Cl), 

ieJ(x) i£l{x) 

where the last three relationships hold since x G f2j for each i € /(a;). This justifies inclusion (|3.11l) . 
In the case of I(x) = 0, we arrive at inclusion (|3.13[) by the first row of the above relationships and 
the normal cone definition (|2.7I) . 

Assertion (ii) is justified similarly to the proof of Thcorcm l3.2f ii') by using the results of assertion 
(i) and the well-known fact that the projection operator for a closed and convex set in a Hilbert 
space is single-valued. A 

Observe that in Theorem 13. 3[ in contrast to Theorem 13. 2[ we do not impose the well-posedncss 
requirement. In fact, under the assumptions of Thcorcm l3.3f ii') it holds automatically; see Corol- 
lary 1.106]. Note also that in finite-dimensional spaces X we always have the Frechet differentiability 
of the distance function at out-of-set points with unique projections (see, e.g., [16[ Exercise 8.53]), 
and so we can deal in the proof of Theorem 13. 3f i) directly with the cost function in the generalized 
Heron problem (|1.1[) . without considering the auxiliary problem (|3.14[) . However, in Hilbert spaces 
this approach requires additional and unavoidable assumptions on the projection continuity; see [5J 
Corollary 3.5]. In finite dimensions the projection continuity and Frechet differentiability of the 
distance functions actually follows from the projection uniqueness, while it is not the case in Hilbert 
spaces as shown in [5l Example 5.2]. Observe to this end that neither uniqueness nor continuity of 
projections is required in Theorem 13.31 

On the other hand, the next result shows that for the unconstrained version of (|1.1[) . i.e., for the 
generalized Fermat-Torricelli problem [13] with disjoint sets Hi, the projection nonemptiness at a 
local optimal solution automatically implies the projection uniqueness in arbitrary Hilbert spaces. 
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Proposition 3.4 (projection uniqueness at optimal solutions). Let x be a local optimal 
solution to problem (|1.4[) in a Hilbert space X with Q = X and n" =1 f2i = 0. Assume that x Qi as 
i = 1, . . . , n. Then the fulfillment of the condition Tl(x; VL{) ^ for all i = 1, . . . ,n implies that the 
projection set H{x; fli) is a singleton whenever i G {1, . . . , n}. 

Proof. Since I(x) = for the first index set in (|3.f 0|) . it follows from the proof of Theorem I3.3f i) 
with f2 = X that for every uji G 11(5;; Qi) as i = 1, . . . , n we have the equality 

n _ 

Picking any Qj, say Or, let us check that the set n(x; fii) is singleton. Indeed, take two projections 
G n(x; fii) and fix arbitrary projections uii G n(5; $7^) for i = 2, . . . , n. Then from (1 3 . 1 5 [) 
we get the relationships 

__ n__ __ n__ 

X — ^1.1 X — LOi X — Wl,2 



d(x;fii) ^d(x;fti) d(x;f2i) ^d(x;fij)' 

which imply that u)i,i = and thus complete the proof of the proposition. A 

Observe that if x belongs to one of the sets 0, as i = 1, . . . , n, the conclusion of Proposition 13.41 
docs not generally hold even in finite dimensions as it is demonstrated by the following example. 

Example 3.5 (nonuniqueness of projections at solution points). Let X = JR 2 in the setting 
of Proposition [331 let be the unit circle of M 2 , and let fi 2 = {(0,0)}. Then x = {(0,0)} is a 
solution of the Fermat-Torricelli problem generated by fli and fi 2 , but the projection H(x;£li) is 
the whole unit circle. It is also clear that any point inside of the unit circle other than (0, 0) is also 
a solution to this problem, and II(x; fli) is a singleton for both i = 1,2, which is consistent with the 
result of Proposition ^. 41 



The observation made in Proposition [331 allows us to improve the optimality conditions obtained 
in p~3l Corollary 4.1] for the generalized Fermat-Torricelli problem. 

Corollary 3.6 (improved optimality conditions for the generalized Fermat-Torricelli 
problem with three nonconvex sets in Hilbert spaces). Let n = 3 in the framework of 
Theorem YS. 31 where Oi , O2 , and Q3 are pairwisely disjoint subsets of X and = X. The following 
alternative holds for a local optimal solution x € X with the sets Ai(x) defined by (|3.12p : 

(i) The point x belongs to one of the sets fij, say f2i. Then for any at G Ai(x) as i = 2,3 we 
have the relationships 

(02,03) < — 1/2 and — ai — as G N(x; Qi). 

(ii) The point x does not belong to all the three sets f2i, Q2, and O3. Then A;(x) = {a^} for all 
i = 1,2,3 and we have 

(ai,a,j) = -1/2 for i^j as i,j e {l,2,3}. 

Conversely, suppose that the sets fli, i = 1,2,3, are convex and that x satisfies either (i) or (ii). 
Then it is a global optimal solution to the problem under consideration. 

Proof. In case (i) for any ai G Aj(x) as i = 2, 3 take G II(x; fi^) such that 

X — LOj 

■ = 2,3. 



d(x; fit) ' 

Since x G fii, we have the relationships 

3 3 



\x - lu 1 \\ + \\x - w 2 || — ^ d(x; Qj) < &i) < d(x; Oi) + \\x - w 2 \\ + \\x - lu 3 \ 



i=l i=l 
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whenever x is near x. Thus a; is a local optimal solution to the problem 

minimize q(x) := d(x; fix) + \\x — a>2 1| + \\x — t*>3 jj . (3.16) 

Employing the generalized Fermat rule in (|3.16p and then the aforementioned sum rule for regular 
subgradients gives us by using the well-known formula for the regular subdifferential of the distance 
function (see, e.g., [9| Corollary 1.96]) that 

G dq(x) = dd(x; fix) + a 2 + a 3 = N(x; fix) H B + a 2 + a 3 . 

The latter implies therefore that 

-a 2 -a3 eiV(5;Oi) with 1 1 a-2 + a.3 1 1 < 1 - 

The rest of the proof follows the lines of that in p~3l Corollary 4.1]. Assertion (ii) and the converse 
statement are derived similarly from Proposition l3.4l and the proof of [HI Corollary 4.1] by the same 
procedure, which thus allows us to fully justify the corollary. A 

From now on in this section we concentrate on the distance function version (jl.ip of the gen- 
eralized Heron problem while paying the main attention to deriving efficient forms of optimality 
conditions for (|1.1[) under additional structural assumptions on the constraint set fi. In what follows 
in this section we impose the nonintersection condition 

fi n fi, = for all i = 1, . . . , n (3.17) 

on the sets fi and fii in (jl.ip . which is specific for the (constrained) generalized Heron problem. In 
this case we obviously have I(x) = for the first index set in (|3.10[) whenever x £ fi, and so the sets 
Ai(x) are calculated by 

Mx) = *~ 1 J li n? i \ i = l,...,n, (3.18) 
a(a:;fij) 

in the Hilbert space setting under consideration. 

To proceed, for any nonzero vectors u, v £ X define the quantity 

(u, v) 

co8(u,v) ■- 



' IN -IN 

and, given a linear subspace L of X, recall that 

L 1 - := {x* 6 X\ (x*,v) = for all v 6 L). 

We say that ft C X has a tangent space L = L(x) at x if L 1 - = N(x; Q). Note that for any affinc 
subspace ft C X parallel to a linear subspace L the tangent space to O at every x £ fi is L. 

Next we derive verifiable necessary and sufficient conditions for optimal solutions to (jl.ip in 
Hilbert spaces provided that the constraint set admits a tangent space at the reference point. 

Proposition 3.7 (optimality conditions for the case of constraint sets with tangent 
spaces). Consider the generalized Heron problem (jl.ip under condition (|3.17p in Hilbert spaces. 
The following assertions hold: 

(i) Let x £ Q be a local optimal solution to (jl.ip . let Ai(x) be computed in (|3 . 18[) where 

n(x; fit) 7^ for i = 1, . . . ,n, and let Q admit a tangent space L{x) at x. Then for any ai(x) £ Ai(x), 
one has 

n 

2^ cos ( a i(^)i v ) — for every v £ L(x) \ {0}. (3.19) 

i=i 

(ii) Let all the sets fi^, i = 1, . . . ,n, be convex. Then Ai(x) = {ai(x)} and condition p.!9p with 
the tangent space L(x) for Q is necessary and sufficient for the global optimality of x, in (jl.ll) . 
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Proof. To justify (i), observe by the assumptions made and the definition of the tangent space L(x) 
to VL at x that 

N(x;n) = L x = {v G X\ (v,x) = for all x <E L(x)}. 



By Theorem 13.31 for any a,i(x) G Ai (x), one has 

n 

e^2a i {x)+L ± (x), 



which implies in turn that 



\^Qi(x),vj = for all v <E 



i=i 



Since x ^ Hi by (|3.17p . we have due to (|3.18[) that ||aj(a:)|| = 1 for i = 1, . . . , n, and hence 

Y n (a / ( ^'? n = whenever « G L(x) \ {0}. 

Thus we arrive at the the necessary optimality condition (|3.19l) . 

To justify (ii), observe that the implication "=>." follows directly from assertion (i) of the 
theorem, since the sets Ai(x) are singletons for i = 1, . . . , n in this case. The oppositive implication 
"<=" follows from Theorem 13. 3f ii) by taking into account the special structure of the normal cone 
N(x; ft) = L J -(x). This completes the proof of the proposition. A 

We have the following specification of optimality conditions in Proposition 13 . 71 when the tangent 
space therein is finitely generated. 

Corollary 3.8 (optimality conditions for the case of finitely generated tangent spaces). 

Let L(x) = span{«i, . . . , v s } with Vj ^0 as j — 1, . . . , s in the setting of Proposition 13.71 Then 
condition (|3.19l) in all of its conclusions is equivalent to 

n 

cos(aj, Vj) = for all j = l,...,s. (3.20) 

i=l 

Proof. We obviously have that (|3.19[) => (|3.20[) . To justify the converse implication, set a := Y27=i °* 
and observe by Vj ^ as j = 1, . . . , s and ||a ? || = 1 as i = 1, . . . , n that ()3.20j) yields (a, Vj) = for 
all j = 1, ... , s. Picking further an arbitrary vector v G L(x) \ {0}, we arrive at the representation 

s 

v = X i v i 
i=i 

with some Xj G M. It gives by linearity that (a,v) = Y^j=i ^j( a i v j) — 0j which yields (|3 . 19f) and 
completes the proof of the proposition. A 

The next result concerns the generalized Heron problem for two nonconvex sets in Hilbert spaces 
with a one-dimensional structure of the regular normal cone to the constraint. 

Proposition 3.9 (necessary conditions for the generalized Heron problem with two non- 
convex sets in Hilbert spaces). Consider problem (jl.ll) for two sets (n = 2) in Hilbert spaces 
under the noninter 'section condition (|3.17[) . Let x <E Q be a local optimal solution to (jl.ll) such that 
N(x; Q) = spanjw} with some u ^ and that H(x; f2j) ^ for i = 1,2. Then for any ai(x) G Ai(x) 
as i = 1,2 we have the conditions: 

either a\{x) + 02(2;) — or cos (ai(x), v) = cos (a,2(x), v). (3-21) 
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Proof. It follows from Theorem 13. 3f i) in this setting that 

- oi(x) - o 2 (x) £ N(x; Q) for any a^x) G A { (x), i = 1,2. (3.22) 

Denoting for simplicity a, := Oj(x) as i = 1, 2 and taking into account the assumed structure of the 
regular normal cone to H, we get that (|3.22p is equivalent to the following: 

either a\ + a 2 = or a\ + a 2 — Xv with some A =/= 0. 

Let us show that the latter condition implies that cos(ai,v) = cos(a 2 ,v). Indeed, in this case we 
have 1 1 di|| = j|ai|| = 1, which gives by the Euclidean norm on X that 

A 2 |M| 2 = IK + a 2 \\ 2 = IKf + ||a 2 || 2 + 2(oi,a 2 > = 2 + 2(a 1; a 2 ). 

This implies in turn the relationships 

(oi, Xv) = (Xv — a 2 , Xv) 

= X 2 \\v\\ 2 -X{a 2 ,v) 

= 2 + 2(01,02) - A(a 2 , u) 

= 2(a 2 , a 2 ) + 2(ai, a 2 ) — X(a 2 ,v) 

= 2(a 2 + ai,a 2 ) — X(a 2 ,v) 

= 2(Xv, a 2 ) — X(a 2 ,v) = (a 2 , Xv), 

which yield that (ai,w) = (a 2 ,v) since A ^ 0. By taking into account that ||ai|| = ||a 2 || = 1 and 
v ^ 0, we conclude that cos(eti, u) = cos(a 2 , v) and thus complete the proof. A 

Observe that sufficient optimality conditions in the form of Proposition 13.91 do not hold even 
in convex settings. The next result provides slightly modified conditions, which are sufficient for 
optimality in the case of the convex generalized Heron problem on the plane. 

Proposition 3.10 (characterizing optimal solutions for the generalized Heron problem 
with two convex sets). Let the sets fii and Q 2 be convex in the setting of Proposition^^ and 
let a t := ttj(S) as i = 1, 2. Then the modification 

either a\ + a 2 = or \a\ ^ a 2 and cos(ai, v) = cos(a 2 , v)\ , (3.23) 

of the necessary condition (|3.2ip is sufficient for the global optimality of x G f2 in (| 1 . 1 1) when 
X = R 2 . 

Proof. To justify the sufficiency of conditions (|3.23[) for the optimality of x in (jl.ip . we need to 
show — by taking into account Theorem l3.3f ii') and the assumed structure of the regular normal cone 
to fl that the relationships in p.23[) imply the fulfillment of 

- 01 - 02 6 N(x; 0) = span{u}. (3.24) 

When — ai — a 2 = 0, inclusion (|3.24[) is obviously satisfied. Consider the alternative in (|3.23[) when 
Oi 7^ a 2 and cos(ai,w) = cos(a 2 ,t>). Since we are in Ft 2 , represent ai = (xi,yi), a 2 = (x 2 ,y 2 ), and 
v = (x,y) with two real coordinates. Then the equality cos(ai,i>) = cos(a 2 ,t;) can be written as 

xix + yiy = x 2 x + y 2 y, i.e., (xi — x 2 )x = (y 2 - yi)y- (3.25) 

Since o ^ 0, assume without loss of generality that y 0. By the equivalence 

IKH 2 = \\a 2 \\ 2 ^ xj + y 2 = x 2 +y 2 

we have the equality (x± — x 2 ){x\ + x 2 ) = (y 2 — yi)(y 2 + yi), which implies by (|3.25[) that 

y(xi - x 2 ){xi + x 2 ) = x(xi - x 2 ){y 2 + yi). (3.26) 
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Note that x\ ^ x 2 , since otherwise we have from (|3.25[) that yi =1/2, which contradicts the condition 
&i 7^ o-2 in (|3.23l) . Dividing both sides of (|3.26[) by x\ — X2, we get 



y(x t + x 2 ) = x{y 2 + yi), 

which implies in turn that 

y(ai + a 2 ) = y(x 1 +x 2 ,yi+ y 2 ) = (x{yx + y 2 ),y(yi + y 2 )) = (yi + y 2 )v. 
In this way we arrive at the representation 

2/1 + 2/2 
a-i + a 2 = v 

y 

showing that inclusion p.24[) is satisfied. This ensures the optimality of x in (jl.ip and thus completes 
the proof of the proposition. A 

We conclude this section by a simple example showing how the results obtained allow us to 
completely solve a direct generalization of the classical Heron problem in M 2 , where the constraint 
straight line is replaced by a convex set. 

Example 3.11 (complete solution of a a convex set extension of the Heron problem on 
the plane). Consider problem where f2 is the epigraph of the nonsmooth convex function 

y = \x\ in M 2 , and where Q\ and are two points (xi,yi) and (x 2 ,y 2 ) that do not lie on Q. This 
problem admits optimal solutions due to Proposition I3.1f ii). To solve it, we are going to employ 
appropriate necessary optimality conditions obtained above. Observe first that the normal cone to 
f2 at (0, 0) is given by 

AT((0,0);O) = {(x,y) e 1R 2 \ y< -\x\} 

while the classical normals at other points of are calculated trivially. Using this, we can easily 
check that if the points (xi,yi) and (x 2 ,y 2 ) belong to the region 

{(x,y) eM 2 \ y < 

then the origin x = (0, 0) is the only point that satisfies the necessary optimality condition from 
Theorem 13. 3f i) written now as: 

( x j y ^ ) 

— a\ — a 2 G N(x; Q) with = — as i = 1, 2. 

If the points {xi,yi) and (22,2/2) belong to another region 

{(x,y) eM 2 \x> \y\}, 

then the problem also has a unique optimal solution constructed by connecting the reflection point 
of (xi, yi) through the line y = x and {x 2 , y 2 ). 

4 Subgradient Algorithm in the Generalized Heron Problem 

In this section we develop a subgradient algorithm for the numerical solution of the generalized Heron 
problem (jl.4l) for finitely many convex sets and convex constraints in the finite-dimensional Euclidean 
space M rn . These are our standing assumptions for the rest of the paper. Recall that n(a;; f2) denotes 
the (unique) Euclidean projection of x to f2 while H^.(x) stands for the generalized/minimal time 
projection (|3.1j) of this point to the target sets ili in (|1.4|) . Here is the algorithm whose various 
implementations are presented in the next section. 
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Theorem 4.1 (subgradient algorithm for the generalized Heron problem). Let S ^ be 

the set of optimal solutions to problem (|1.4[) . Picking a sequence {ak}keiN of positive numbers and 
a starting point x\ G Q, consider the algorithm 

n 

Xk+i = Ilffffc - Qifc qjk] Oj , k = 1,2,..., (4.1) 

wii/i an arbitrary choice of vectors 

q ik G -dp F {oj ik - x k ) ("1 iV(w ife ; Q;) /or some G (x fe ) i/ x k £ (4.2) 
via the Minkowski gauge (|3.2[) and with q^ := otherwise. Assume that 

OO OO 

a k = oo and £ 2 := a 2 , < oo. (4-3) 
fe=i fe=i 

TTien i/ie iterative sequence {x k } in (|4.1[) converges to an optimal solution of problem (|1.4[) and i/ie 
numerical value sequence 

V k :=mm{T(x j )\j = l,...,k} (4.4) 
converges to the optimal value V in this problem. Furthermore, we have the estimate 

- dOr i; S) 2 + L 2 £LX 



v k -v< 



2 



where < L < oo is a Lipschitz constant of the function T(-) /rom (|1.4p on i? m . 

Proof. We know that the value function T(-) in f| 1 .4f) is convex and globally Lipschitzian on M m . 
Employing |12[ Theorems 7.1 and 7.3], the convex subdiffcrcntial of the minimal time functions (|1.3[) 
at any point x k is computed by 

C N(x k ;Qi) n{veX\ (T F (-v) < 1} if x k G ^, 
OT&(a*) = ^ (4.5) 
[ iV(w ifc ; n<) n [ - dp F (bJ ik - x k )] if £ fe ^ Qj, 

where u;^ G HQ.(x k ) is an arbitrary generalized projection vector for i G {l,...,n} and fc G iV. 
Recalling now the subgradient algorithm for minimizing the convex function T(-) in (|1.4[) subject to 
x G f2, we construct the iteration sequence by 



a; fe+ i = U[x k - a/c^/c; fi) with v k G dT(x k ), k = 1,2, (4.6) 

It follows from the convex subdiffcrcntial sum rule of Theorem I2.1f ii) that 



Vk 



^2 qik with G dT^.{x k ) 



i=l 



for the subgradicnts in (|4.6|) . Substituting the latter into (|4.6p gives us algorithm (|4.1[) with 
Qjfe satisfying (|4.2I) . Then all the conclusions of the theorem are derived from the so-called "square 
summable but not summable case" of the subgradient method for constrained convex functions 
under the conditions in (|4.3|) : see [IJE] for more details. A 

In the case of F = IB, the closed unit ball in M m , we are able to provide a more explicit 
algorithm to solve the distance function version (|1.1[) of the generalized Heron problem with now 
uniquely defined vectors qik in (|4.1I) . 
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Corollary 4.2 (explicit subgradient algorithm for the distance version of the general- 
ized Heron problem) . Consider the distance function specification of the generalized Heron 
problem under the assumptions of Theorem 14.11 Then all the conclusions of this theorem hold with 
qik in (|4.1|) calculated by 

{0 if x k e Cli, 

x k -TL(x k ;fk) , n (4J) 

— T — if x k f fti. 

Proof. As follows from the proof of Theorem 13.31 in the case of problem (jl.ll) the vectors qik from 
(|4.2|) are uniquely determined and reduce to (|4.7[) . A 

The next corollary specifics algorithm (|4.1[) in the case of balls for the distance function version 
(|l.lj) of the generalized Heron problem. 

Corollary 4.3 (subgradient algorithm in the case of multidimensional balls). Consider 
problem with Oj = B(ci, r^) C R™ as i = 1, . . . ,n. Then the quantities qik in Theorem W.W are 
uniquely calculated by 

{0 if \\x k -Ci\\ < n, 

x k - Cl „ ( 4 - 8 ) 

and the corresponding values Vk are evaluated by formula (|4.4[) with 



n 

T( Xj )= Yl (ll^-cill -»•*)■ ( 4 - 9 ) 



Proof. Formula (|4.8p directly follows from (|4.7p due to the projection representation 



n(x k ]£ii 



\\Xk-Ci\\ 



in the case under consideration. It is easy to see furthermore that the value function in (|4.4|) reduces 
to (Oil in this case. A 



5 Implementation of the Subgradient Algorithm 



The final section of the paper is devoted to implementations of the subgradient algorithm from 
Theorem 14.11 and its specifications to solve the generalized Heron problem in a number of underly- 
ing examples of their own interest. Let us start with a two-dimensional problem involving a ball 
constraint in the setting of Corollary [ 




-10 -8 -6 -4 -2 



MATLAB RESULTS 


k 


Xk 


v k 


1 


(-1,4) 


44.58483 


10 


(-1.07737,3.61433) 


44.36969 


100 


(-1.07779,3.61332) 


44.36969 


1000 


(-1.07779,3.61331) 


44.36969 


10,000 


(-1.07779,3.61331) 


44.36969 



Figure 1: A Generalized Heron Problem for Balls with a Ball Constraint. 
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Example 5.1 (two-dimensional Heron problem for balls with ball constraints). Consider 
the generalized Heron problem (|l.ip for balls in M 2 subject to a given ball constraint. Let Cj = (di,bi) 
and r, as i = 1, ...,n be the centers and the radii of the balls under consideration, and let 
c = (xo,yo) and r be the center and radius for the given ball constraint fl. The subgradient 
algorithm is given by (|4.1[) . where the projection P(x 7 y) :~ Tl((x, y); f2) is computed by 



p (x, y) = (v x +x ,v y + y ) with v x 



r(x - x ) 



r(y - yo) 



(x - x ) 2 + (y - yo) 



2 ' 



x /(x-x ) 2 + (2/-yo) 2 ' 



and where the quantities qik and 14 are calculated in Corollary [ 

To specify the calculations, take the ball constraint fl with center (—2,4) and radius 1. The 
sets Qi, i = 1, . . . , 6, are the balls with centers (-10, 0), (-1, 8), (2, -4), (7, 6), (7, 1), and (8, -3) 
and with the same radius r = 1. The MATLAB calculations performed by algorithm (|4.1[) with the 
sequence ctk = 1/fc satisfying (|4.3[) and the starting point x\ — (—1,4) are presented in Figure 1. 
Observe that the numerical results indicate points on the ball constraint with the optimal solution 
x w (-1.07779,3.61331) and the optimal value V « 44.36969. 

The next example concerns the generalized Heron problem with square constraints. 




MATLAB RESULTS 


k 




v k 


1 


(-1,-4) 


41.23881 


50 


(0.89884,-3) 


37.32496 


100 


(0.95169,-3) 


37.32091 


150 


(0.97352,-3) 


37.31974 


200 


(0.98595,-3) 


37.31920 


250 


(0.99413,-3) 


37.31890 


300 


(1.00000,-3) 


37.31872 


350 


(1.00000,-3) 


37.31872 



Figure 2: A Generalized Heron Problem for Balls with a Square Constraint. 



Example 5.2 (generalized Heron problem with square constraints). Consider the imple- 
mentation of algorithm (|4.1j) for problem (|1.1|) using a MATLAB program with the square constraint 
fl of center (a, b) = (0, —4) and short radius r = 1 and with the balls Qi as i = 1, . . . , 6 centered 
at (-7,-3), (0,5), (-4,0), (2,-4), (6,0), and (6,7) with the same radius 0.5. Note that the projection 
P(x, y) — n((x, y); CI) is calculated by 

a + r, b + r) if x — a > r, y — b > r, 

x,b + r) if \x — a\ < r, y — b > r, 

a — r, b + r) if a; — a < — r, y — b > r, 

a ~ r i J/) if i — a < — r, y — b| < r, 

P(x, j/) = < (a — r, b — r) if x — a < — r, y — b < — r, 

x, 6 — r) if |x — a| < r, y — b < — r, 

a + r, b — r) if x — a>r, y — b< —r, 

a + r,b) if x — a > r, |y — b\ < r, 

x,y) if (x,y) G Q. 
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The quantities a^. and Vk arc given by Corollary |4.3l In Figure 2 we present the results of calculations 
performed by the subgradicnt algorithm (|4.1I) for the sequence ak = 1/k and the starting point 
x\ = (—1,-4). Observe that the computed optimal solution is x ~ (1.00000,-3.00000) and the 
optimal value is V w 37.31872. 

Prior to the calculations in two next examples concerning the generalized Heron problem 
for squares in M 2 we formulate a specification of Theorem 14.11 in a general setting of such a type. 
Recall that a square in M 2 is of right position if the sides of this square arc parallel to the x-axis 
and the y-axis, respectively. 

Corollary 5.3 (subgradient algorithm for the generalized Heron problem squares tar- 
gets). Consider problem fll.lj) in M 2 , where each target set fli is a square of right position with 
center Cj = (<2j, 6j) and short radius r.j as i — l,...,n, and where the constraint Q is an arbi- 
trary closed and convex set. Denote the vertices of the i th square by Vu = (a; + r,,6j + 7"j), V2i = 
(at - n,bi + n), V3i = (a, - n,bi - r»), vu = (a, + n,bi - r»), and Zei a; fc = (in, a; 2 fe). TTiera i/ie 
quantities in Theorem 14.11 are computed by 

if \xik - ai| < n and |x 2 fc - 6i| < 7\, 
if xife — aii > r, and X2fc — 6, > n, 



9ifc = < 





- Vu 


II** 


-Vli\\ 


a;* 


- Vu 


\\Xk 


-V 2i \\ 


Xk 


- V3i 


\\Xk 


-«3i|| 


Xk 


— VU 



\\Xk - V4i 

(0,1) 
(0,-1) 
(1,0) 
(-1,0) 



if xik — at < —n and x 2k — h > n, 
if xik —cu< —ri and x 2k — bi < —n, 

if xik — ai > ri and x 2 k — bi < —ri, 

if \xik - Oi| < n and x 2 k - &i > n, 
if \xik - Oi| < Ti and x 2k - b t < -r % , 
if xik — ai > ri and \x 2 k — fo\ < n, 



if Xik — ai < —n and \x 2 k — b»| < n 
for all i = 1, . . . ,n and k E N with the corresponding quantities Vk defined by (|4.4p . 



Proof. This statement follows from Corollary 14.21 by a direct calculation of the projection from an 
out-of-set point to each square fJi in formula (|4.7I) . A 

Now we present the results of MATLAB calculations in the case of straight line constraints in 
the setting of Corollary [ 



Example 5.4 (generalized Heron problem for squares with line constraints). Consider 
the generalized Heron problem (|1.1[) for squares of right position in M 2 subject to a straight line 
constraint tt. Let Cj = (a^, bi) and r, as i = 1, . . . , n be the centers and short radius of the squares 
fli under consideration. Denote by Vu = (eij + fj, bi + r,), «2i = (a, — r.;, 6j + r^), = (a, — r^, 6j — 
r i), u 4i = ( a i + r i, &t — r i) l ne vertices of the i th square, and let v = [s, h] and p = (xo, 2/o), be the 
direction and point vectors of the given line f2. Then the projection P(x,y) = H((x,y);H,) in the 
the subgradient algorithm (|4.1[) is calculated by 



P(x, y) = (xq + st, y + ht) and t 



s(x - x ) + /?.(?/ - j/p) 

S 2 + /l 2 



while the quantities qik and for all i = 1, . . . , n and k G N are given by Corollary 15.3 
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MATLAB RESULTS 


k 




v k 


1 


(-1,6) 


42.8838 


100 


(-1.0826,6) 


42.8821 


1000 


(-1.0896,6) 


42.8821 


100,000 


(-1.0938,6) 


42.8821 


1,000,000 


(-1.0944,6) 


42.8821 


5,000,000 


(-1.0946,6) 


42.8821 


10,000,000 


(-1.0946,6) 


42.8821 



Figure 3: A Generalized Heron Problem for Squares with a Line Constraint. 



In Figure 3 we present the results of calculations by algorithm (|4.ip with a k = 1/k and the 
starting point x\ = (—1,6) for the case above with the line constraint defined by v = [1,0] and 
p = (1,6) and the squares Oj as i = 1, . . . , 5 centered at (—6, —9), (—5, 4), (0, —7), (1, 0), and (8, 8) 
with the same short radius r=l. Observe that the calculated optimal solution is x ~ (—1.0946,6) 
and the optimal value is V « 42.8821. 

The next example concerns the generalized Heron problem (jl.ip for squares in right position 
with a ball constraint on the plane. 




MATLAB RESULTS 


k 


Xk 


Vk 


1 


(5,-2) 


54.41891 


10 


(3.51379,-1.33835) 


53.05740 


100 


(3.41230,-1.21623) 


53.04403 


1000 


(3.39607,-1.19475) 


53.04364 


100,000 


(3.39279,-1.19033) 


53.04363 


600,000 


(3.39271,-1.19022) 


53.04363 


1,000,000 


(3.39271,-1.19021) 


53.04363 


1,200,000 


(3.39270,-1.19021) 


53.04363 


1,400,000 


(3.39270,-1.19021) 


53.04363 



Figure 4: A Generalized Heron Problem for Squares with a Ball Constraint. 



Example 5.5 (generalized Heron problem for squares with ball constraints). By taking 
into account the previous formulas for algorithm (|4.ip . we provide the following calculations con- 
cerning the generalized Heron problem (|1.1[) with the ball constraint fi centered at (5, 0) and radius 
2 and the squares £1;, i = 1, . . . , 8 of right position with the centers (—2, 4), (—1, —8), (0, 0), (0, 6), 
(5, —6), (8, —8), (8, 9), and (9, —5) and the same short radius r = 0.5. Figure 4 presents the results 
of calculations for algorithm (|4.1|) with the sequence a k = 1/fc and the starting point x\ = (5, —2). 
Observe that the obtained numerical results give us the optimal solution x « (3.39270,-1.19021) 
and the optimal value V « 53.04363. 

Now let us illustrate applications of the sub-gradient algorithm from Theorem 14. II to solving the 
generalized Heron problem (|1.4[) formulated via the minimal time function with dynamics sets F 
different from the ball. First we consider the dynamics F described by the closed unit diamond 

F := {(x u x 2 ) € M 2 \ |xi| + \x 2 \ < 1}. (5.10) 
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In this case the corresponding Minkowski gauge (|3.2I) is given by the formula 

Pf(xi,x 2 ) = \x\\ + \x 2 \. 



(5.11) 



The following proposition provides an explicit calculation of a sub-gradient of the minimal time 
function (|1.3[) generated by the diamond dynamics (|5.10p and a square target in M 2 . We further use 
this calculation in implementing algorithm (|4.1[) with the corresponding selection of qik in (|4.2j) . 



Proposition 5.6 (subgradients of the minimal time function with diamond dynamics). 

Let F be the closed unit diamond in 1R 2 , and let f2 be the square of right position centered at c = (a, b) 
with short radius r > 0. Then we can calculate a subgradient v(xx,x 2 ) G dT^ (xi, x 2 ) of the minimal 
time function Tq at (xx,x 2 ) ^ by 



v(x ly x 2 ) = < 



(1,0) 


if \x2 — b\ < r 


xi > a + r 


(-1,0) 


if \x2 — b\ < r 


x\ < a — r 


(0,1) 


if \xi — a\ < r 


X2 > b + r 


(0,-1) 


if \xi — a\ < r 


X2 < b — r 


(1,1) 


if x~i > a + r, 


X2 > b + r, 


(-1,1) 


if xi < a — r, 


X2 > b + r, 


(-1,-1) 


if xi < a — r, 


X2 < b — r, 


(1,-1) 


if Xi > a + r, 


X2 < b — r, 





if (xi,X2) e SI 





(5.12) 



Proof. By [12l Theorem 7.3] we have the relationship 

dT£ (x) = N(ui; 0) n [ - dp F (u> - x)] for any Q G Il£(x) 



(5.13) 



between the subdiffercntials of the minimal time function at x f2 and the corresponding Minkowski 
gauge. In the setting under consideration it is easy to find the minimal time projection Hq(xi,x 2 ) 
of a given vector (xi,x 2 ) G M 2 to the square Q. Furthermore, the convex subdiffcrcntial of (|5,11[) 
at (xi,x 2 ) is computed by 



dpF(xi,X2) = < 



[-1,1] x [-1,1] 


if 


(xi,x 2 ) = (0,0), 


[-1,1] X{1} 


if 


X\ = 0, X2 > 0, 


[-1,1] x{-l} 


if 


Xi = 0, X2 < 0, 


{1}X[-1,1] 


if 


Xl > 0, X2 = 0, 


{-l}x[-l,l] 


if 


X\ < 0, X2 = 0, 


{1}x{1} 


if 


Xl > 0, X2 > 0, 


{I}x{-1} 


if 


Xl > 0, X2 < 0, 


{-1}x{1} 


if 


Xl < 0, x 2 > 0, 


{-I}x{-1} 


if 


Xl < 0, X2 < 0. 
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The rest of the proof is a direct verification that the vector v{x\,X2) from (|5.12[) belongs to the set 
on the right-hand side of (|5.13|) and hence to dT^(x\ : X2)- A 

Proposition 15.61 and the previous considerations lead us to the following realization of the sub- 
gradient algorithm (|4.1[) . 

Corollary 5.7 (subgradient algorithm for finitely many squares and diamond dynamics 
in the generalized Heron problem). Consider problem (|1.4[) generated by the diamond dynamics 
(15.1 Op and n squares fli of right position in IR 2 . Let c, = (aj,6,) and as i = l,...,n be the 
centers and the short radii of the squares under consideration, and let Vu = (a^ + + fj), 

V2i = (a. - n, h + n), v 3l = (oj - n,bi - n), and vu = (a, + r.;, ^ - r<) &e the vertices of the i th 
square. Denoting Xk — (xik,%2k) i n algorithm (|4. 1|) . we compute the quantities qik as follows: 



lik = < 



' 


if \xik - 


- a 


< n and \x 2 k 


- bi\ < n, 


(1,1) 


if x lk - 


a. 


> n and X2k — 


bi > n, 


(-1,1) 


if xik - 


04 


< —rj and X2k 


- bi > n, 


(-1,-1) 


if xik - 


(I, 


< — n and x 2 k 


- bi < —n 


(1,-1) 


if xi k - 


(I, 


> n and x 2 k — 


bi < -n, 


(0,1) 


if \ssik - 


II 


< n and X2k - 


-bi > n, 


(o,-i) 


if \xi k - 


- a 


< n and x 2 k - 


-bi < —n, 


(1,0) 


if xik - 


a, 


> n and \x 2 k - 


bi\ < n, 


- (-1,0) 


if xik - 


a, 


< — n and \x 2 k 


-h\ <n 



for all i G {1, . . . , n} and k G IN . 

Proof. It follows from Proposition l5.6[ comparison between the right-hand side of (|4.2[) and formula 
(|5.13p . and the square calculations of Corollarv l5.3l A 

Now we implement the results of Corollary 15.71 to solve the generalized Heron problem of the 
above type with ball constraints. 




-6 -4 -2 2 4 6 8 10 



MATLAB RESULT 


k 


Xk 


v k 


1 


(1,-2) 


34 


10 


(1.98703,-0.83947) 


32.01297 


100 


(1.99987,-0.98385) 


32.00013 


1,000 


(2.00000,-0.99838) 


32.00000 


10,000 


(2.00000,-0.99984) 


32.00000 


50,000 


(2.00000,-0.99997) 


32.00000 


100,000 


(2.00000,-0.99998) 


32.00000 


150,000 


(2.00000,-0.99999) 


32.00000 


200,000 


(2.00000,-0.99999) 


32.00000 



Figure 5: A Generalized Heron Problem for Squares with a Ball Constraint with Respect to "Sum" 
Distances. 



Example 5.8 (generalized Heron problems with diamond dynamics for squares and ball 
constraints). Consider problem (|1.4j) with the diamond dynamics (|5.10|) for squares 17^ as i — 
1,...,6 of right position in Ft 2 with the centers at (-5,-3), (-4,0), (2,3), (4,-5), (5,6), and (8,-1) 
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and the same short radius 1 subject to the ball constraint f2 centered at (— 1, 1) and radius 1. The 
results of calculations by the subgradicnt algorithm (14.11) with a k = 1/k and the starting point 
x\ = (1,-2) are presented in Figure 5. Observe that the obtained optimal solution is the point 
x « (2.00000, -0.99999) on the ball constraint with the optimal value V « 32.00000. 

The following example is a modification of the previous one for the case of square constraints. 



MATLAB RESULT 


k 




v k 


1 


(-1,2) 


61 


10 


(1,0.57897) 


54.5 


100 


(1,0.48990) 


54.5 


1000 


(1,0.50100) 


54.5 


1500 


(1,0.49933) 


54.5 


2000 


(1,0.49950) 


54.5 


2500 


(1,0.49960) 


54.5 


3000 


(1,0.50000) 


54.5 


3500 


(1,0.50000) 


54.5 



Figure 6: A Generalized Heron Problem for Squares with a Square Constraint with Respect to 
"Sum" Distances. 



Example 5.9 (generalized Heron problems for squares with diamond dynamics and 
square constraints). Consider the generalized Heron problem (|1.4j) with the diamond dynam- 
ics (|5.10|) for the squares fl; G M 2 as i = 1, . . . , 7 of right position centered at (—5, —3), (—9, 1), 
(0, 6), (2, —3), (6, 8), (5, —5), and (9, 1) with the same short radius 1 subject to the square constraint 
f2 of right position centered at (0,1) with the short radius 0.5. The calculations presented in Figure 6 
are performed for the sequence a k = 1/k in (|4.1[) and the starting point X\ = (— 1, 2). The obtained 
optimal solution is the point x sw (1, 0.50000) on the square and the optimal value is V « 54.50000. 

Next we consider the generalized Heron problem (|1.4[) with the square dynamics F = [—1,1] x 
[— 1, 1] on the plane. The corresponding Minkowski gauge is now given by 



Pf{xi,x 2 ) 



First we calculate a subgradient v(xi, x 2 ) € 9Tq x 2 ) of the cost function in (|1.4[) at any (xi, x 2 ), 
which is further used for a specification of algorithm (|4.1I) in this setting. 



Proposition 5.10 (subgradients of minimal time functions with square dynamics and 
square targets). Let F = [—1, 1] x [— 1, 1], and let Q be the square of right position in M 2 centered 
at c = (a,b) with short radius r > 0. Then a subgradient v(xi,x 2 ) £ 9Tj^(xi,X2) of the minimal 
time function Tq at (xi,x 2 ) is computed by 



j(xi,x 2 ) = < 



(1,0) 


if \X2 


-b\ 


< Xl 


- a, 


Xi > a + r 


(-1,0) 


if \X2 


- b\ 


<a- 


Xl, 


xi < a — r 


(0,1) 


if \x\ 


-a 


< X2 


-b, 


X2 > b + r, 


(0,-1) 


if \xi 


-a 


<b- 


X2, 


X2 < b — r, 





if (xi 


x 2 ) 


e ft. 







(5.14) 



Proof. It is given in [13j Proposition 5.1]. A 

As a consequence of the proposition above, we calculate the quantities in algorithm (|4.ip for 
the corresponding version of the generalized Heron problem. 
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Corollary 5.11 (subgradient algorithm for the generalized Heron problem with square 
dynamics). Consider problem (jl.4l) for the square dynamics F = [—1,1] x [—1,1] and the square 
targets fli as i = l,...,n of right position in M 2 . Denote by Ci = (a,-,6j) and r» the centers 
and the short radii of the squares f2j under consideration, and let the vertices of the i th square be 
Vu = (Oj + r-j, bi + n), v 2l = (a l - r,-, + n), v 3l = (a* - r u b l - r,-), and v 4l = (a* + n,bi - r,). 
TTien the quantities qik in algorithm (|4.1[) o/ Theorem 14.11 m t/iis setting along the iterative sequence 
Xk = (iCifc, ^2fe) flre calculated for all i G {1, . . . , n} and k £ IN by 

(1,0) i/ a^2fc - &i| < xik - at and x lk > a* + n, 

(—1,0) if \x 2 k - h\ < <M — Xik and x lk < a t - n, 

(0, 1) if \xik — o»| < X2fc — &, and X2k > h + r», 

(0, -1) i/ lain, - Oi| < 6, - a- 2 fe and a; 2 fe <h-n, 

(0, 0) otherwise. 



qik 



and formula 

A 



Proof. It follows from Proposition l5 . 1 ll comparison between the right-hand side of 
dSHH) , and the square calculations of Corollarv l5.3l 

The following two examples present implementations of the subgradient algorithm realization 
from Corollary 15.111 in the generalized Heron problem under consideration with square and ball 



MATLAB RESULT 


k 


Xk 


v k 


1 


(-4,3) 


26.25000 


10 


(-3.12500,1.04603) 


24.37500 


100 


(-2.99136,1.00070) 


24.25068 


1,000 


(-3.00133,1.00133) 


24.25000 


10,000 


(-2.99996,1.00001) 


24.25000 


15,000 


(-3.00013,1.00007) 


24.25000 


20,000 


(-3.00000,1.00000) 


24.25000 


25,000 


(-3.00001,1.00001) 


24.25000 


30,000 


(-3.00001,1.00001) 


24.25000 



Figure 7: A Generalized Heron Problem for Squares with a Ball Constraint with Respect to "Max" 
Distances. 




MATLAB RESULT 


k 




v k 


1 


(5,0) 


35 


10 


(4.00062,0.03519) 


33 


100 


(4.00000,0.00038) 


33 


200 


(4.00000,0.00010) 


33 


400 


(4.00000,0.00002) 


33 


600 


(4.00000,0.00001) 


33 


800 


(4.00000,0.00001) 


33 


1,000 


(4.00000,0.00000) 


33 


1,200 


(4.00000,0.00000) 


33 



Figure 8: A Generalized Heron Problem for Squares with Ball Constraint with Respect to "Max" 
Distances. 
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Example 5.12 (generalized Heron problem with square dynamics, targets, and con- 
straints). Consider the implementation of the algorithm from Corollarv l5.11l in problem (|1.4p with 
the square constraint f2 of center (-3,2) and short radius 1 and the target square sets f2j as z = 1, . . . , 5 
of centers (-8,6), (-6,-2), (-1,8), (-1,-7), and (2,6) with the same short radius r = 0.75. In Figure 7 
we present the results of calculations by (|4.1|) with at = 1/fc and the starting point x\ = (—4,3). 
The optimal solution here is x k, (-3.00001, 1.00001) and the optimal value is V « 24.25000. 

Example 5.13 (generalized Heron problem with square dynamics and targets and with 
ball constraints). Consider the implementation of the subgradient algorithm from Corollary 15.111 
in problem (|1.4[) with the square dynamics, the square targets fli as i = 1, . . . , 6 of centers (-5,-8), 
(-4,5), (0,0), (8,7), (9,3), and (7,-3) with the same short radius r = 0.5, and with the ball constraint 
f2 of center (5,0) and radius 1. The presented calculations are performed by (|4.1j) with oik = l/k and 
the starting point x\ = (5,0); see Figure 8. The obtained optimal solution is x (4.00000,0.00000) 
with the optimal value V « 33.00000. 

Our last example concerns a three-dimensional distance version of the generalized Heron problem 
(|l.lj) for cubes of right position in M 3 subject to a ball constraint. 




MATLAB RESULT 


k 


Xk 


v k 


1 


(5,0.5,-6) 


51.58786 


10 


(4.23949,1.52680,-4.79680) 


47.19028 


100 


(4.23948,1.53023,-4.79546) 


47.19026 


1,000 


(4.23948,1.53024,-4.79546) 


47.19026 


10,000 


(4.23948,1.53024,-4.79546) 


47.19026 


100,000 


(4.23948,1.53024,-4.79546) 


47.19026 


1,000,000 


(4.23948,1.53024,-4.79546) 


47.19026 



Figure 9: A Generalized Heron Problem for Cubes with Ball Constraint in Three Dimensions. 



Example 5.14 (generalized Heron problem for cubes with ball constraints). Consider 
problem (jl.ip for cubes Oj as i = 1,...,6 of right position in M 3 with the centers (8,-4,3), 
(—2, —6, 3), (3, —2, 2), (—4, —5, —6), (—3, 1, 1), and (3, 7, —5) and the same short radius 1 subject to 
the ball constraint fl of center (5,2,-6) and radius 1.5. The projection P((x, y, z); fl) and quantities 
qik in algorithm (|4.1j) are calculated similarly to Example [5Tj Figure 9 presents the implementation 
of the subgradient algorithm (|4.1[) with ak = 1/k and the starting point x\ = (5, 5, —6). As we see, 
the optimal solution calculated here up to five significant digits is x. « (4.23948, 1.53024, —4.79546) 
and the optimal value is V « 47.19026. 

We conclude the paper by the following three observations. 

Remark 5.15 (extensions and other location problems). 

(i) Note that the approach and results of this paper can be easily extended to the weighted 
version of the generalized Heron problem (|1.4[) : 

n 

minimize T(x) := V^/ijTj^(a;), subject to x S ft, (5.15) 

i=l 

where /i,- > as i = 1, . . . , n are given weights. Since we have 

d( fH Tl)(x)= IM dTg i (x) 
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for both convex and nonconvex subdiffercntials used in this paper, it is straightforward to derive 
counterparts of the qualitative and numerical results obtained above for the case of the weighted 
generalized Heron problem (|5.15|) . For example, the equation 



replaces the one in (|3.19[) for all the corresponding results. 

(ii) Our variational approach can be used to solve a variety of other facility location problems. 
In particular, the following smallest intersecting ball problem can be naturally formulated and in- 
vestigated by using the above tools of variational analysis and generalized differentiation: given n 
nonempty closed subsets fli C X, i = 1, . . . ,n, find a point lona given set f2 and the smallest 
number r > such that the ball with center at x and radius r has nonempty intersection with all 
the sets fij as i = 1, . . . , n. This problem is modeled as follows: 



We intend to address this and other facility location problems in our future research. 

(iii) For some results in the Hilbert space setting of Section 3, it is possible to use the prox- 
imal normal cone instead of the Frechet normal cone. However, we use the Frechet normal cone 
consistently for the simplicity of presentation. 
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